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Elliptic flow in transport theory and hydrodynamics 
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We present a new direct simulation Monte-Carlo method for solving the relativistic Boltzmann 
equation. We solve numerically the 2-dimensional Boltzmann equation using this new algorithm. 
We find that elliptic flow from this transport calculation smoothly converges towards the value 
from ideal hydrodynamics as the number of collisions per particle increases, as expected on general 
theoretical grounds, but in contrast with previous transport calculations. 
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Ultrarelativistic nucleus-nucleus collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC) have been argued to 
create a "perfect liquid" , with an extremely low viscos- 
ity The essential piece of evidence is the large mag- 
nitude of elliptic flow at RHIC which is as large as 
predicted by ideal-fluid models (which assume zero vis- 
cosity). Elliptic flow is an azimuthal asymmetry in the 
momentum distribution of particles, projected onto the 
plane transverse to the beam direction {z axis): in a colli- 
sion between two nuclei with non-zero impact parameter, 
more particles are emitted parallel to impact parameter 
{x axis) than perpendicular to it (y axis). This asymme- 
try results from the almond shape of the overlap region 
between the colliding nuclei, which is transformed into 
a momentum asymmetry by pressure gradients Q: mi- 
croscopically, elliptic flow results from the interactions 
between the produced particles, and is therefore a key 
observable of the dense matter produced at RHIC. 

In the microscopic language of particle physics, ideal 
fluid and low viscosity translate into small mean free 
path of a particle between two collisions or, equivalently, 
large rescattering cross sections between the "partons" 
created in a collision. The description of the system in 
terms of partons is itself questionable at the early, dense 
stage of the collision, but it is nevertheless a helpful, in- 
tuitive picture. The natural question which arises then 
is: how large must the partonic cross section be in order 
to achieve ideal-fluid behavior, i.e., local thermal equili- 
birum? How many partonic collisions are needed? 

In this paper, we address this issue by solving numer- 
ically a relativistic Boltzmann equation, and comparing 
the results with relativistic hydrodynamics. It is well 
known that the Boltzmann equation reduces to hydro- 
dynamics when the mean free path is small (see '3| for a 
rigorous proof in the relativistic case). Numerically, how- 
ever, it has been found [sj that hydrodynamics produces 
larger elliptic flow (by 30-40%) than the Boltzmann equa- 
tion. In this paper, we address this issue using a differ- 
ent method. A possible explanation for the discrepancy 
found in Q is suggested at the end of this paper. 

The primary limitation of the Boltzmann equation is 
that it only applies to a dilute system, where the mean 
free path of a particle is much larger than the distance 



between particles, so that one need only consider two- 
body collisions, many-body collisions occurring at a much 
lower rate. There is no reason to believe that the RHIC 
liquid is dilute: interactions are nonperturbative, so that 
both the mean free path and the distance between par- 
ticles are of order 1/T, where T is the temperature. 
Clearly, the Boltzmann equation cannot be used to di- 
rectly simulate a heavy-ion collision; nevertheless, it has 
the potential of giving us a grasp on deviations to ther- 
malization. 

The relativistic formula for the collision rate between 
two beams of particle densities rii and n^, and arbitrary 
velocities Vi and V2 can be written as 



coU 



dtd^x 



cmin2 V|vi - V2p - |vi x V2P/c2, 



(1) 



where a is the scattering cross section. There are several 
methods for implementing this collision rate in Monte- 
Carlo algorithms. The most widely used method 0, Q 
is the ZPC algorithm [^]: this algorithm treats particles 
as hard spheres, which collide when their distance in the 
center-of-mass frame is smaller than r = y/ a/ir. An 
alternative method is the stochastic collision algorithm 
proposed by Xu and Greiner [l3| • 

Here, we implement a new algorithm, where relativis- 
tic effects are incorporated in a physically transparent 
way. The second term under the square root in Eq. ([1]) is 
the relativistic correction. It is simply understood, in the 
case of colliding hard spheres, as a geometrical effect re- 
sulting from the Lorentz contraction of the spheres. We 
take this contraction into account, in the Monte-Carlo 
simulation, by replacing spheres with oblate spheroids, 
whose polar axis is the direction of motion. As with the 
other algorithms, Lorentz covariance and locality are bro- 
ken in the sense that collisions occur instantaneously at a 
finite distance. However, these violations are small if the 
system is dilute @, which is required by the Boltzmann 
equation. The difference with the ZPC algorithm is that 
the collision time is determined directly in the laboratory 
frame, not in the center-of-mass frame. 

For sake of simplicity, we consider a two-dimensional 
gas of massless particles in the transverse plane {x,y). 
We thus neglect an important feature of the dynamics 
of heavy-ion collisions, the fast longitudinal expansion. 
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FIG. 1: Picture of a collision between a particle 1 of size 
o"2d ~ '2r and a pointlike particle 2. The impact parameter d, 
as defined in Eq. (|3]), is negative. 



As will be shown below, this is not crucial for elliptic 
flow, which is essentially unaffected by the longitudinal 
expansion. Initial conditions for the N particles are gen- 
erated randomly according to a gaussian distribution in 
coordinate space, and an isotropic, thermal distribution 
in momentum space: 
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t-y > Rx and Pt = \Jp'i: + Py- The temperature 
r(x, y) is determined as a function of the local particle 
density per unit area, n2d according to n2d oc T^, the 
equation of state of a massless ideal gas in 2 dimensions. 

Particles interact via 2 — > 2 elastic collisions. We as- 
sume for simplicity that the total elastic cross section 
a2d is independent of the center-of-mass energy s. (In 
2 dimensions, the cross section has the dimension of a 
length.) This is implemented in a Monte-Carlo calcula- 
tion by treating each particle as a rod (due to Lorentz 
contraction) of length r — \(J2d perpendicular to the ve- 
locity. Equivalently, for each pair of colliding particles, 
one can assume that particle 1 has length 2r and particle 
2 is pointlike (see Fig. [1]). The impact parameter of the 
collision is 



d=- 



Gz ■ ((X2 - Xi) X (V2 - Vl)/c) 
1 — Vl • V2 /c^ 



(3) 



This quantity is Lorentz invariant. Our collision algo- 
rithm is deterministic: the scattering angle in the center- 
of-mass frame 9* is determined as a function of d, as 
in classical mechanics. The results presented here are 
obtained with an isotropic differential cross section, i.e., 
9* = 7r(l — d/r). This prescription ensures that particles 
move away from each other after the collision and do not 
collide again. 

We now define two dimensionless numbers relevant to 
this problem. The average particle density per unit sur- 
face n2d is [13 
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AirR^Ry 



(4) 
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The typical distance between two particles is n^^ , while 
The mean free path of a particle between two collisions 
is A = \/{(T2dn2d)- The ratio of these two lengths is the 



dilution parameter D: 
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(5) 



As explained above, applicability of Boltzmann theory 
requires D -^1. In addition, locality and covariance re- 
quire that the interaction length be much smaller than 
the mean free path [1, i|. In two dimensions, the in- 
teraction length is a2d, and CT2d/A = D^: locality and 
covariance are thus recovered in the limit Z? <C 1. 

The second dimensionless number is the Knudsen num- 
ber Kn, which characterizes the degree of equilibration by 
comparing A with the system size R. The latter quantity 
can be measured by any average of Rx and Ry . A natural 
choice for elliptic flow is [ll[ 



-1/2 



The Knudsen number is then defined as 
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(6) 



(7) 



Hydrodynamics is the limit Kn 1, while the limit 
Kn ^ 1 corresponds to free-streaming particles. The 
values of D and Kn can be tuned by varying the cross 
section t72d and the number of particles N . 

Note that the mean free path is a local quantity, which 
depends on space-time coordinates: in particular, it in- 
creases as the system expands. Strictly speaking, D and 
Kn defined by Eqs. ([5]) and ([7]) are initial values. Dimen- 
sional analysis suggests that they are the relevant control 
parameters for this problem. 

Solving the Boltzmann equation in the hydrodynamic 
limit requires both D <C 1 and Kn ^ 1. This in turn 
requires a huge number of particles N in the Monte-Carlo 
simulation: inverting the above equations, one obtains 
(assuming Rx — Ry for simplicity) 



N 



Stt 



D2Kn^ 



(8) 



Since the computing time scales grows with N like iV3/2 
one cannot implement arbitrarily small values of D and 
Kn. Instead, we choose to study numerically the depen- 
dence of elliptic flow on D and Kn, and extrapolate to 
the hydro limit D = Kn = 0. 

Fig. [2] presents our results for the average elliptic flow 
V2 = (cos2(/)) = {(pI -pI)/{pI +P^)). The aspect ratio 
of the initial distribution is Ry = 1.5Rx, and the Monte- 
Carlo simulation has been pushed to very large times 
(t = lOOi?), so that all collisions are taken into account. 
Quite naturally, V2 decreases monotonically to as Kn 
increases towards the free-streaming limit Kn — > -|-oo. 
For a given value of D, the variation with Kn is very well 
described by the simple formula proposed in fllj |: 



V2 



1 + Kn/Kno' 



(9) 
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FIG. 2: Variation of the elliptic flow «2 with the Knudsen 
number, Kn, for several values of the dilution parameter D. 
The statistical error on each point is 5v2 — 7 x 10"*. For each 
value of D, Monte-Carlo results are fitted using Eq. (|9}. 
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FIG. 3: Time dependence of the average elliptic flow V2 from 
the transport model, from the corresponding 2-dimensional 
hydro calculation, and from usual 3-dimensional hydro with 
Bjorken longitudinal expansion. 



where the parameters V2 and Kng are fit to our Monte- 
Carlo results. V2 is the limiting value of V2 when Kn — > 0, 
expected to coincide with V2 from hydrodynamics in the 
limit D —> 0. Surprisingly, the value of V2 depends very 
little on D: v!^ = 0.102 ± 0.003. This means that the 
hydrodynamic limit is more general than the Boltzmann 
equation and applies even if the system is not dilute. Un- 
like V2, the parameter Kuq strongly depends on D. We 
do not have a simple explanation for this dependence. 
However, only the limit D <^ 1 has a well-defined phys- 
ical interpretation, as it corresponds to the Boltzmann 
equation. For larger values of D, locality and causality 
are broken, and the physical interpretation of the results 
is less clear. For D 1, our fit gives Kuq = 0.70 ± 0.03. 

An independent hydro calculation with the same initial 
conditions was done using the same code as in Ref. 
For sake of consistency 5], the equation of state of the 
fluid is that of a two-dimensional ideal gas (i.e., the equa- 
tion of state of a dilute g modeled by the Boltzmann 
equation) , whose velocity of sound is Cg = c/ |13i| . The 
average V2 at t ^ lOOR is V2^'^'° = 0.101 ± 0.003 (the er- 
ror bar in the hydro calculation is due to the fact that 
the calculation is pushed to very large times). This is 
compatible with the value V2 obtained from the Boltz- 
mann calculation: Boltzmann transport theory and ideal 
hydrodynamics agree in the limit Kn — > 0, as expected. 

The Knudsen number is closely related to the average 
number of collisions per particle ficou. From the defini- 
tion, Eq. ([7|), one expects that the number of collisions 
per particle is ^ 1/Kn. The proportionality constant 
can be computed exactly for a gaussian distribution in 
the limit Kn ^ 1 [l3| , which yields the relation 

2A^coii _^V2_ f 2e \ 1 1.6 

rirnn = = ; K Kn ~ , 

N ^3/2^1^77 \l + ej Kn' 

(10) 

where the factor 2 means that each collision involves 2 



particles, e = {Ry — R'^)/{Iiy + R'^) is the initial eccentric- 
ity, and K(x) is the complete elliptic integral of the first 
kind. We used this formula to check numerically that 
our algorithm produces the right number of collisions for 
large Kn. In addition, our numerical results show that 
the product ncouKn is remarkably constant for all values 
of Kn, so that Eq. (fTO| holds within a few percent. Us- 
ing Eq. ([9|), this shows that an average of 2.3 (resp. 9.1) 
collisions per particle are required to achieve 50% (resp. 
80%) of the hydrodynamic "limit" on V2- 

Fig. [3] compares the time-dependence of elliptic flow 
in hydro and in the transport model. The subtle point 
is that the convergence of V2{t) towards hydro as Kn 
is not uniform: for a given value of Kn, deviations 
from hydro are large at early times, and tend to decrease 
afterwards. More precisely, it can be shown, following 
the same methods as in Ref. li*], that the early-time 
behavior is V2 oc in the transport model, and V2 oc in 
hydro. This is compensated by the late-time behavior: V2 
decreases slowly at large times in hydro (not seen in the 
figure), while V2 from the transport model stays constant 
after the last collision has occurred. 

Can our 2-dimensional results be used in the context of 
heavy-ion collisions? The essential difference in 3 dimen- 
sions is the fast longitudinal expansion due to the strong 
Lorentz contraction of the colliding nuclei. Elliptic flow, 
however, is a purely transverse observable which is little 
affected by the longitudinal expansion: Fig. [3] displays a 
comparison between the average elliptic flow computed 
in 2-dimensional hydro and in 3-dimensional hydro with 
Bjorken longitudinal expansion (the initial time in this 
calculation is tq — i?/4, corresponding to tq ~ 0.4 fm/c 
in a semicentral Au- Au collision) [l^ . Both yield similar 
results for the magnitude and the time-dependence of ti2. 

One can reasonably expect that deviations from hydro 
are similar in 2 dimensions and in 3 dimensions. The 
number of collisions per particle, ncoii, is not the right 
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TABLE I: Expected values of V2 for semi-central Au-Au col- 
lisions at RHIC. For each value of the cross section, we quote 
the value of the equivalent Debye-screened, leading order 
QCD cross section ^. 



Isotropic a [mb] 


Debye-screened a [mb] 


Kn 




3 


7 


0.72 


0.49 


8 


19 


0.27 


0.72 


20 


47 


0.11 


0.87 



quantity for carrying out the comparison (neither the 
transport opacity fS]): in 3 dimensions, ncoii is large at 
early times (with a fixed partonic cross section, ricoii di- 
verges like luTg'^ as the initial time tq goes to 0), but 
these collisions do not produce much elliptic flow (see 
Fig. [3]). As argued in Ref. [ll|, the Knudsen number 
should be evaluated at the time when elliptic flow de- 
velops t ^ R/cs, with Cs — l/\/3 in the quark-gluon 
plasma phase. The mean free path is A = I /an with 
n = {l/ct){l/S){dN/dy) and S = iTrR^Ry, hence 

1 IdNcs , , 

dN/dy is the total (charged-l- neutral) multiplicity. For a 
semi-central Au-Au collision at RHIC, {1/ S){dN/dy) ~ 
0.8 mb~^. With this definition of Kn, we expect that 
Eq. ([9]) should hold approximately in 3 dimensions, with 
Knp ~ 0.7. Table [J gathers numerical estimates ob- 
tained using the same values of a as in Ref.Q. For 
CT = 20 mb, corresponding to a QCD cross-section of 
47 mb, the transport result should be only ~ 10 — 15% 
below hydro. 

For this value of the cross section, Molnar and Huovi- 
nen find that V2 from the transport calculation is lower 
by 30% than V2 from hydro Q . Furthermore, the depen- 
dence of their results on cr is at variance with our results. 
Using Eqs. (O and pT|) . the deviations from ideal hydro 
should decrease as 1 /cr as cr increases. This is very gen- 
eral: dissipative effects are expected to be linear in the 
viscosity ry, which scales like l/a. The discrepancy with 
hydro should be at least twice smaller with cr = 47 mb 
than with a = 20 mb, which is clearly not the case in 
Fig. 1 of Ref. In our opinion, this cannot be at- 
tributed to dissipative effects. 

The origin of the problem may be the dilution con- 
dition D ^ 1, which is not satisfied in previous trans- 



port calculations. The trick to reduce D in the ZPC 
cascade algorithm is the "parton subdivision technique" : 
one multiplies the number of particles iV by a large num- 
ber I, and one divides the partonic cross section cr by 
Z, so that the Knudsen number Kn is unchanged. The 
distance between particles in three dimensions is n~^/^, 
which replaces n^j^ in Eq. ([5|) . Taking parton subdi- 
vision into account, one obtains D — an^/^l^'^/^ . The 
average particle density in a semicentral Au-Au collision 
at time t = R/cs is 2.5 fm"^ [ll|. With cr = 20 mb and 
I — 180 0, D ~ 0.65, which is not very small compared 
to unity. The situation is worse at early times due to 
the longitudinal expansion: D ~ 1.1 at tq — 0.6 fm/c. 
For Kn — 0.11, carrying out the calculation with D ~ 1 
leads to overestimate the deviation from hydro by at least 
a factor of 2 (see Fig. [2]). The reason why previous cal- 
culations were done with such large values of D is that 
small values are hard to achieve numerically: in 3 dimen- 
sions, Eq. ([8]) is replaced with N oc D~^Kn~^ , and the 
computing time grows like A^^/"^ oc D^'^Kn^*. This was 
our motivation for studying first the 2-dimensional case. 

In summary, we have implemented a new algorithm 
for solving the relativistic Boltzmann equation. Wc have 
studied the convergence of the relativistic Boltzmann 
equation to relativistic hydrodynamics. Since this re- 
quires both a dilute system (D <C 1) and a small mean 
free path (Kn <C 1), which is costly in terms of computer 
time, we have restricted our study to two dimensions. 
This preliminary study has only addressed the average 
value of the elliptic flow; differential results as a function 
oipt, as well as results on v^, will be presented in a forth- 
coming publication. We have shown that the average 
elliptic flow computed in the transport algorithm con- 
verges smoothly towards the hydrodynamics value as the 
strength of final-state interactions increases. Isotropic 
partonic cross sections of 3 mb and 20 mb should pro- 
duce respectively ~ 50% and ^ 90% of the hydro value 
for V2- 
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